Physics of Proximity Josephson Sensor 
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We study the proximity Josephson sensor (PJS) in both bolometric and calorimetric operation 
and optimize it for different temperature ranges between 25 mK and a few Kelvin. We investigate 
how the radiation power is absorbed in the sensor and find that the irradiated system is typically in 
a weak nonequilibrium state. We show in detail how the proximity of the superconductors affects 
the device response: for example via changes in electron-phonon couphng and out-of-equilibrium 
noise. In addition, we estimate the applicability of graphene as the absorber material. 

PACS numbers; 



I. INTRODUCTION 

As a class of high sensitivity infrared detectors, super- 
conducting sensors operating at cryogenic temperatures 
are increasingly popular. For different design choices of 
the device, superconductivity guarantees greatly reduced 
heat leakage^^ drastic drop of resistance at the critical 
temperature^ or negligible Johnson noise when operat- 
ing in dissipationless regime.— The detector characteriz- 
ing noise equivalent performance (NEP) values are thus 
expected to be as small as 10~^^ W/-\/Hz at T « 1 K 
for the state-of-the-art detectorsi^ Depending on whether 
the sensor operates in continuous or single-photon detec- 
tion mode, it can be used in astrophysical space appli- 
cations, security applications, spectroscopy of ultrafast 
quantum phenomena, optical communications, quantum 
cryptography or fast digital circuit testing. 

There exist several measures to characterize the per- 
formance of a detector. Furthermore, their relative im- 
portance depends on the specific application in question. 
The requirements for a good detector include low noise, 
sufficient speed and good impedance matching between 
the receiver antenna and the detector element. Different 
detectors such as transition edge sensors and kinetic in- 
ductance detectors excel in different areas and operation 
regimes. When considering any one particular high-end 
application, the choice of the detector needs to be deter- 
mined on a case-by-case basis, depending on the qualities 
and properties required. 

Very recently, a model for a mesoscopic kinetic- 
inductance detector has been proposed.^ This device 
is based on the coherent superconducting proximity ef- 
fect, in contrast to some other superconductor-normal- 
metal-superconductor (SNS) detectors, ''-li where the su- 
perconductors are only used to block the energy outdiffu- 
sion from the radiation absorber and transport between 
the superconductors remains incoherent. The Proximity 
Josephson Sensor (PJS) under consideration consists of 
a long SNS junction (length larger than the supercon- 
ducting coherence length) , where the temperature in the 
normal island changes under irradiation. This change in 
temperature causes an exponential change in the Joseph- 
son critical current between the two superconductors, 
which, in turn, greatly modifies the kinetic inductance 



FIG. 1: (Color onhne) Illustration of the detector with normal 
(N) and superconducting (S) components. Incident radiation 
is coupled to the normal-conducting absorber via the antennas 
in contact to the superconductors. 



of the sensor. This can be read with high accuracy us- 
ing for example SQUID or resonant circuit based readout 
methods. As a result, performance values such as NEP 
- 7 X 10-2" W/%/Hi (T = 200 mK) and signal-to-noise 
ratio ^ 10^ (/ ^ THz) are expected, when the detec- 
tor is operated in continuous and single-photon modes, 
respectively. In addition to favorable performance ex- 
pectations, the advantages of this type of a detector in- 
clude versatility, i.e., the fact that the device can be easily 
fabricated to optimize detection at the chosen operating 
temperature. Moreover, the response can be measured 
from a well-defined and spatially small region where tem- 
perature changes. The smallness of this region is relevant 
since the sensitivity depends on electron-phonon interac- 
tion, the strength of which is proportional to the volume. 
Finally, it is possible to get a good understanding of the 
physics of the device, as we show below. 

Here we detail the proposition made in Ref. 0, for- 
mulate a model to explain the process of radiation cou- 
pling to the sensor and also take a closer look on how the 
proximity effect induced in the normal metal by the su- 
perconductors affects the performance of the detector. A 
number of distinguishable effects are due to the proxim- 
ity induced minigap in the quasiparticle excitation spec- 
trum of the normal island. Some examples of these are 
the reduced electron-phonon scattering and the reduced 
electron heat capacity. While considering the underly- 
ing physics of the detector, we also take on a practical 
approach to determine the required operating conditions 
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more closely than has been done in Ref. The analysis 
contains some aspects general to all detectors with a sim- 
ilar hot-electron mechanism, especially when we study 
the coupling of the incident radiation to the absorber and 
the mechanism by which the radiation power is converted 
to a thermal signal. When the operating conditions are 
varied, different physical realizations in terms of detector 
materials may become favourable and therefore we also 
consider whether graphene could be used as the radiation 
absorber. The potential benefits in this case stem from 
the low volume and the subsequent decrease in intrinsic 
thermal noise of the two-dimensional sheet of carbon. 



II. DETECTOR 

1. Basic operating principle 

The detector consists of a diffusive metal wire (N) con- 
nected to two superconductors (S) through transparent 
contacts, as portrayed in Fig.[T] Due to the proximity ef- 
fect, superconducting correlations penetrate the normal 
wire and a steady supercurrent can be sustained between 
the superconducting terminals, in the absence of an ex- 
ternal voltage. The maximum value of this supercurrent, 
Ic, is strongly dependent on the electron temperature Tg 
in the wire in the temperature range i?Th < ksTe A, 
where A is the energy gap in S and Etyi = HD/P is 
the Thouless energy of the N wire with diffusion coeffi- 
cient D. In short, incoming radiation increases electron 
temperature in the absorber, this increase in temperature 
makes the Josephson critical current Ic{Te) decrease, and 
the decrease in Ic translates to a change in the kinetic 
inductance of the system, Lk ~ l//c- The change in Lk 
can then be measured by one of the standard techniques 
discussed below. 

Since the detector responds to incident radiation with 
elevation of temperature, it is essential that the temper- 
ature range where the response is sufficiently large to be 
measured is within reasonable limits. In the diffusive 
wire, this range is determined by the energy minigap 
Eg{(f)) induced in the wire due to the proximity effect. 
This minigap depends on the phase difference (j) between 
the superconductors and for <^ = it obtains its max- 
imum Eg(0) = 3.1i?Th-^ When the wire is long enough 
so that the overlap of S order parameters in the mid- 
point is small, i?Th ^ A and the operating regime of 
the detector is determined to be well below the critical 
temperature of the superconductors. On the other hand, 
strong superconductors are needed to prevent outdiffu- 
sion of excited electrons into the leads before they relax. 
This outdiffusion grows larger with larger i?Th which is 
proportional to the diffusion constant D. For small D, 
the critical current becomes lower and consequently the 
operating temperature decreases. The exact regime de- 
pends on the measurement accuracy for the supercurrent 
and the desired frequency range of operation. 



2. Detailed discussion 

The remainder of this section is devoted to a more de- 
tailed description of the operating procedure. Here and 
below, we strive for using simple analytical expressions 
for the essential physical properties of the detector. In- 
homogeneous superconductivity in diffusive conductors 
can be described with the Usadel equationi^ To describe 
the proximity induced effects, we solve the Usadel equa- 
tion numerically and calculate the response coefficients 
such as heat conductance between electrons and phonons 
from the solutions J^iii This is required, for example, to 
determine the modified density of states for the absorber 
material. We then write down and use approximate ex- 
pressions drawn from these exact numerical results. Sec- 
ondly, we use a model where the effect of the electromag- 
netic radiation coupling to the absorber is assumed to 
show as an increase in the electron temperature of the 
absorber. The question of radiation coupling to super- 
conducting elements is a long-standing problem which 
we unfortunately cannot answer in detail in the course of 
this discussion. Nevertheless, we expect that the result 
of this coupling can be described by a hot-electron model 
in our case, where the frequency of the incident photons 
v ~ THz is much larger than the frequency correspond- 
ing to the edge of the energy minigap in the absorber 
Eg/h ~ 5 GHzJ^ We detail our reasoning in Sec. Ill Al 

The critical Josephson current for a diffusive wire can 
be determined from the Usadel equation. For E^h ^ 
ksTe < A, it ifM 

' ~ (3 + 2^/2)ei^Jv V ^Th 1^ V ^Th / 

GAnkBTe I { ^ \ 
~ (3 + 2V2)eRN ^^"^ V C^/ ■ 



Here, we have defined £,n = \/ hD/2TTkBTe, and the 
normal state resistance of the wire, Rn ~ P^/A — 
l/Npe'^DA, is expressed in terms of the wire resistiv- 
ity p, wire cross section A, and normal-state density of 
states Np at the Fermi level, in addition to the quantities 
already defined. The numerical solution for the critical 
current as a function of electron temperature together 
with the analytical approximation of Eq. ([T|) are shown 
in Fig. [5] for some typical values of material parameters 
to be defined below. The regime of strongest response 
can be seen to come about at temperatures of a few i^xh 
and while Ic saturates at low temperatures, we see that 
values of ksTc > -Exh are already sufficient for a sizable 
response dIc/dTc- This sets the lower limit for opera- 
tion. The upper limit is set by the minimum measurable 
critical current, and here we assume this to be of the or- 
der of a few nanoamps corresponding to some tens of nH 
in the inductive readout j^. In Fig. [2l critical current of 
3 nA is obtained at temperatures of roughly 400 niK, cor- 
responding to 20i<^Th- It is notable that if the linearity 
of response is the preferred characteristic, the detector 
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FIG. 2: (Color online) Critical Josephson current Ic as a func- 
tion of electron temperature Te. We use ratio A/^Jxh = 830 
for E'r-h/kB — 19 mK, corresponding to Nb energy gap 
A = 16 KxkB. 



should be operated at a narrow temperature range of 
STf, ~ _ETh within the chosen operating point. 

Changes in Ic are readily reflected in the kinetic in- 
ductance for a Josephson weak link, which is defined as 
Lk = n/2e{dls/d(j)y^ « V2e(cos(/))-i. The current- 
phase relation for the supercurrent is typically approx- 
imately sinusoidal, Is = IcSm{(j>). The most feasible 
value for to be used in operation of the sensor needs 
to be determined in the experiments and it depends on 
the measurement accuracy for the supercurrent. The 
readout can be performed by using a SQUID or prob- 
ing the resonant frequency of an LC resonator where the 
SNS junction provides the inductance. Specifically, in 
the previous case the magnetic flux noise can be as small 
as 10"'' ^o/Vilz (<I>o is the magnetic flux quantum)^^ so 
that the contribution from the readout can become 1-2 
orders of magnitude smaller than the contribution from 
intrinsic sources and therefore negligible.— The sensitiv- 
ity of the detector is then ultimately determined by the 
intrinsic sources, i.e., energy fluctuations in the absorber. 

The response to the incoming radiation and the cor- 
responding raise in temperature is mostly dependent on 
the volume of the element absorbing the radiation en- 
ergy and the degree of energy relaxation in the electron 
subsystem. We do not exclude the possibility of ther- 
mal transport into the leads but it is preferable to design 
the detector in a way to prevent the excess energy in the 
electron system escaping the island. Then, the electron- 
phonon coupling becomes the dominant mechanism for 
energy relaxation and this is in the normal state charac- 
terized by the electron-phonon coupling constant E. The 
size of the detector is essential in three ways. First, as 
shown below, the intrinsic noise of the sensor scales with 



respect to the volume of the radiation absorbing element 
and secondly, a small absorber responds more strongly 
to small incoming quanta of energy in terms of elevated 
temperature. Moreover, the normal metal bridge has to 
be short enough to be able to support measurable super- 
current. Thirdly, however, the detector cannot be too 
short or the outdiffusion rate of energy grows consider- 
ably. This means that typical scales lie near one micron 
for detector length and range from tens to hundreds of 
nanometers for cross-sectional dimensions. In addition 
to the absorber and the superconducting leads coupled 
to it, a separate antenna is required to couple the in- 
coming radiation to the actual SNS sensor. The input 
impedance of the N wire can be varied to match the an- 
tenna impedance resulting in a good quantum efficiency. 
Most notably, a typical resistance of the N bridge is close 
to 50 il which usually gives the best matching. 



A. Coupling to radiation 

The radiation detector studied in this work is based 
on a small electronic system absorbing the radiation. In 
such systems even a weak perturbation may give rise to a 
notable change in the state of the system. As shown be- 
low, in some cases this state cannot entirely be described 
via a simple increase of temperature. In the following, we 
aim to answer the three questions: (i) how is the energy 
distributed inside the absorber after the initial excita- 
tion, (ii) to which extent can this energy absorption be 
described by an increase in the electron temperature and 
(iii) how does the energy finally escape the system. A 
proper theory of the effects of radiation in a coherent 
SNS system for the frequencies of the order of i?Th is 
quite complicatedi^. However, for frequencies hi^ > 2A 
we expect any coherent effects to be smallji^ Therefore, 
at frequencies close to terahertz it is justified to assume 
that coupling to the proximity device is similar to cou- 
pling to a normal metal. 

To address the question of energy distribution and es- 
cape, we formulate a model based on the Boltzmann 
equation in a steady state, 

= -^rad(e, + /c-ph(e) + /c-c(e) + 

(2) 

— (/(e)-/eq(e, Tbath))e(|e| - A), 
td 

where the electron energy distribution function at en- 
ergy e, /(e), is modified from its equilibrium Fermi-Dirac 
form /Ar(e, Tbath) due to incident radiation of frequency 
u. The collision integrals /rad(e, i^), ^c-ph(e), ^c-c(e) de- 
scribe the interaction between the electron having energy 
e and the radiation field, absorber phonons and other 
electrons, respectively?i^ The radiation collision integral 
is composed of all the processes corresponding to absorp- 
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tion and emission of a photon with frequency 



Iiad — 



^(/(e + M(l - fie)) + fie M(l - fie)) 
- - fie + hu)) - /(e)(l - /(e - hu))) 

-Popt 



{hv) 



:(/(e + H + /(e-M-2/(e)), 



and the incident optical power is Popt- In the absence 
of coherent corrections the prefactor Po^tHhv)^ is inde- 
pendent of energy. For |e| > A, the possibihty of the 
quasiparticle escaping to the superconducting leads is 
described by the diffusion term with the diffusion time 
td = ETh/fi and Heaviside step function 0(|e| — A). 
The latter accounts for the lack of heat diffusion into the 
superconductors for energies below A. The origin of re- 
laxation time for diffusion, td/8, is easily understood by 
noting that the quadratic length dependence of td gives a 
factor 4 for diffusion length L/2 from the midpoint of the 
wire and since there are two possible directions for dif- 
fusion, the escape time is reduced from td by the factor 
8. In the superconductor at |e| > A, the quasiparticles 
are assumed in equilibrium at the bath temperature with 
distribution fsiThnth)- The total power corresponding to 
each of these processes is obtained by multiplying by eNp 
and integrating over energy. 

We solve Eq. ^ to observe that the response of the 
system to incident radiation is almost thermal at typical 
materials parameters and a wide range of input powers 
Popt and frequencies i^. This means that the electron 
distribution function /(e) of the irradiated system can, 
to a good accuracy, be approximated by an equilibrium 
Fermi function with a certain electron temperature Te 
determined by input power and some other parameters. 
Consequently, the supercurrent response is independent 
of the fine details of the nonequilibrium. However, Tg 
cannot be obtained by a simple quasiequilibrium heat 
balance between the radiation power and the electron- 
phonon coupling. In most cases, this means that some 
fraction of the optical power flows out from the system 
in a way that it does not contribute to equilibrium heat- 
ing of electrons and Te < T^, where is the resulting 
temperature in the case where an equilibrium system is 
heated with power Popt and the energy flows fully into 
the phonons, as assumed for example in Ref. 0- There 
are two distinct reasons for this loss of power. Below, we 
detail these processes. 

The different contributions to electron energy balance 
defined in Eq. ^ are shown for our example system 
in Fig. 131 Since electron-electron interaction only mixes 
electrons with different energies and does not remove en- 
ergy from the system, the corresponding power is zero. 
Due to the conservation of energy, the power correspond- 
ing to the radiation contribution then has to equal the 
integral over the sum of the electron-phonon integral and 
the diffusion part weighed by eNp- The relative mag- 
nitude of these two parts gives the fraction of power 
which either escapes to the leads or to the phonons via 



■Radiation contribution 
-Electron-plionon contr. 
■Electron-electron contr. 
■Diffusion contribution 




-1 



FIG. 3: (Color online) Relative contributions to the electron 
energy balance of the system defined in Eq. ([2|. At all en- 
ergies e, the magnitude of the input (radiation contribution) 
must equal to the sum of the magnitudes of the other three 
processes which contribute to energy relaxation. The result is 
obtained for a typical set of operating and materials parame- 
ters and, most notably, hu — 4A and ksTe « 0.09A. 



the electron system in the absorber while only the lat- 
ter results in a detector response by the means of an 
increased temperature. This is first of the two causes 
for reduced thermal power absorption and it is the most 
prominent for frequencies hv > A. When /ij/ 3> A, 
most of the radiation power is distributed at high ener- 
gies where strong electron-phonon relaxation dominates 
outdiffusion from the island and the relative fraction of 
escaping power becomes small. This may be quanti- 
fied within the the hot-electron assumption and relax- 
ation time approximation for electron-phonon interac- 
tion, when le-ph = (/(e, Te) - /Ar(e, rbath))/Te-ph(e). 

As determined numerically, electron-phonon interaction 
is the dominant factor for inelastic relaxation and we may 
simplify Eq. ^ considerably. Since hiy ^ kBT^, /rad is 
almost constant, /rad — sign(e)Popt/(/iJ^)^0(/iz^ — |e|), for 
energies below the radiation frequency. The correspond- 
ing power is divided between outdiffusion and inelastic 
(electron-phonon) relaxation with ratio 



-^c — ph 



Ptot 2C{8T^'e{\e\ - A) +re_ph(e)-i)ede 



for the power absorbed in the detector. The 
energy-dependent electron-phonon time Tc_ph(e) = 
72C(5)/c|7VF/(Se3)iiia, where ({x) is the Riemann zeta 
function, (^(5) « 1.0369, and thus 



ri = 1 



lU0Ci5)k%NF ihv) 



(hi.) 



(3) 
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The ratio approaches unity for hv ^ l\ and has its min- 
imum at hv « 1.3A, independent of other parameters. 

The second cause of reduced thermal absorption is 
due to the relative weakness of electron-electron inter- 
action. It turns out that for a typical small diffusive wire 
the electron-electron interaction becomes negligible com- 
pared to the electron-phonon interaction at energies that 
are generally smaller than the operating scales we are 
interested in, namely A. This is in contrast to some ear- 
lier assumptions about the source of relaxation in sim- 
ilar hot-electron bolometers.^ The absence of electron- 
electron collisions changes notably the relaxation and 
distribution of the incident radiation power as electron- 
phonon collisions are required to balance the radiation 
power especially at high energies. In Fig. |3] this is seen 
as a deviation of the electron-phonon contribution from 
its quasiequilibrium form where a sharp peak at low ener- 
gies decays with width determined by the electron tem- 
perature. Instead, in the radiation-induced nonequilib- 
rium state electron-phonon interaction is strong also at 
energies e ^ T^s- The resulting increase in the elec- 
tron temperature is smaller than expected from simple 
heat balance. We attribute this to the fact that electron- 
phonon relaxation at high energies is much faster than 
at low energies, resulting in a large heat conductance for 
the nonequilibrium component excited by the radiation 
power. Numerically, we observe that the fraction of inci- 
dent power converted to thermal response (in the absence 
of diffusion) r2 ^ l/{hi>) at high frequencies. Compared 
to the frequency dependence, the depence of r2 on the 
radiation power is relatively weak. 
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independent. The results for the case where both the 
outdiffusion (ri) and the distribution of energy through 
inelastic relaxation r2 are taken into account are shown 
in Fig. 2] for our example system. For small hi', the fre- 
quency dependence is complicated due to the competi- 
tion between the two effects but when hiy ^ A, the loss 
of power follows the 1 / i/-law. 
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FIG. 5: (Color online) Nonequilibrium distribution resulting 
from irradiation and a set of relevant equilibrium distribu- 
tions. Notice that the curves for the actual nonequilibrium 
distribution /(e) and the equilibrium distribution f{e,Tcs) 
used to approximate this are nearly on top of each other. 
The difference between these two is magnified in the inset. 
The figure is obtained for the same set of parameters used 
in Fig. [3l In this case the input power Popt is large so that 
the bath temperature Tbath is small compared to the final Te. 
Here r « 0.08, so Te is also clearly smaller than the tem- 
perature of an equivalent equilibrium system heated with the 
same power, . 

In Fig. [5l the nonequilibrium electron distribution is 
compared to related equilibrium distributions. The es- 
sential difference between the nonequilibrium distribu- 
tion and its quasiequilibrium approximation is the slight 
increase in the fraction of electrons at high energies. 
While this changes essentially the electron-phonon col- 
lision integral in Eq. © to balance the high-energy ra- 
diation excitations even for hv ^ A, it has no relevant 
effect on the supercurrent response of the system. We 
therefore use the quasiequilibrium assumption with elec- 
tron temperature Tf, below as a good approximation. 



FIG. 4: (Color online) Fraction r{i/, P) of incident radiation 
power converted into thermal response in a PJS detector. The 
results are given as a function of v for P = 10 fW and as a 
function of P for hu — 2 A (inset). The curves are obtained 
for a typical set of materials parameters defined in Sec. IIIBI 

The effects described by ri and r2 are not completely 



B. Parameters and conditions 

To quantify the discussion, we assume leads of su- 
perconducting Nb and a wire of normal-conducting Ag 
(also Cu is a viable alternative). As discussed in the 
previous section, a large superconducting energy gap A 
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is needed to prevent energy outdifFusion from tlie ab- 
sorber and to guarantee good quantum efficiency. Tlie 
clioice of superconductor is made solely with this in mind 
and it does not affect the performance in any other way 
since it is the minigap Eg instead of A which determines 
the properties of the detector, as long as the condition 
A ^ ksTe, i?Th is satisfied. The lead parameters are 
then Tc~9K,A = 1.4 meV, and for the wire, we assume 
length I = 1 /im, width w — 100 nm, thickness t = 30 nm, 
Np = 10'*'^ J"^ m-3, D = 25 cm^s"^ This gives 
n = IQ-^^ m^, Rj^ « 52n, ~ l.e^icV^ 19 mKxfcs 
and the proximity minigap in the wire at (j) — is Eg fa 
60 mK= 1.2 GHzx/i. The hmit A > E^xh is weU satis- 
fied with A/ETh — 830. Furthermore, we use electron- 
phonon coupling coefficient S = 5 x 10^ Wm^"^ char- 
acteristic for Ag.^'^ We use these values in the examples 
throughout the paper until Sec. IIII C[ where we study 
different detector setups. 



While the acknowledged advantage in hot electron mi- 
crobolometers is the small size of the absorbing element, 
resulting in increased responsivity, the important pa- 
rameter in this particular design is the Thouless energy, 
i?Th — hD/l"^, determined by the diffusion constant D 
and the length I of the absorber. We note that it is 
precisely the large value of i?Th which differentiates PJS 
from the previous SNS-type detectors. For a detector 
of 1 /im in length, the coherent transport between su- 
perconductors persist up to temperatures comparable to 
the minigap Eg ^ 3.1i?Th ^ 60 mK and as discussed 
above, exploiting the decay of this coherence is the in- 
tegral part of the detector operation. Note that the op- 
erating range for the detector extends well beyond Eg 
and can be determined from Fig. [2l In the SNS detec- 
tors of Refs. H and 0, the absorbers are roughly 5 times 
longer and, since Eg ~ i?Th ~ they reside in the 
incoherent regime for temperatures above 10 mK. Con- 
sequently, the sole purpose of the superconductors there 
is to prevent energy outdiffusion from the absorber. On 
the other hand, to observe strong proximity effect the 
absorber film needs to be sufficiently clean as well. For 
ultrathin silver samples with thickness of 50 nm, diffu- 
sivities of 200 cm^/s have been measuredj-^ but to min- 
imize the noise of the detector, it is preferable also to 
minimize the thickness, even in the expense of diffusiv- 
ity increased by the surface scattering. Our objective to 
strengthen the superconducting correlations by requiring 
clean samples is in contrast to the conventional approach 
in hot-electron detectorsii However, in contrast to Ref. d 
where D = 100 cm^ is used, we need to compromise in 
choosing weaker diffusivity to prevent absorbed energy 
from escaping to the superconductors. This corresponds 
to reduction of i^xh by a factor of 4 which consequently 
dulls the proximity effect, decreaces the critical current 
and shifts the operating regime of the detector to lower 
temperatures. 



III. OPERATION 

The device can be designed to operate in two modes de- 
pending on the rate of the incident photons: continuous 
and single-photon detection. These are referred to as the 
bolometric and calorimetric modes, respectively. Below, 
we study the system with particular focus on the effects 
that stem from the superconducting correlations induced 
in the absorber by the proximity effect. For comparison, 
we refer to the case where these correlations are taken 
into account as the coherent case and the one, where the 
correlations are neglected, as the incoherent case. 

A. Bolometer 

When the detector is subject to a continuous radiation 
with power Popt, its temperature is determined by the 
balance equation for incoming and outgoing powers, 

r{i')Popt ^ Pt = Qc-ph- 

Factor r(z/) accounts for the effects related to the cou- 
pling of radiation, discussed in Sec. Ill Bl and for frequen- 
cies hi' > ksTe it is almost independent of temperature. 
The factor r(i/) may also contain effects such as possible 
impedance mismatch in calibration of the antenna circuit 
which determines the quantum efficiency with which the 
dissipation of the electromagnetic energy can be focused 
on the absorbing element N. After the incoming radia- 
tion energy is first absorbed in the electron system, the 
superconducting gap in the leads prevents the heat from 
escaping the island. The system is consequently relaxed 
by the means of energy flow from the electrons to the 
phonons acting as a heat bath. In the normal state, the 
heat fiux is given by— 

Qf_ph = - T,\,h), (4) 

and the proximity effect can be described by using an 
approximate expression^! 

Qc-ph = Qo-ph^ ^ (5) 

The exponential decay of the interaction is character- 
ized by a constant, T*, and it is due to the presence 
of the minigap Eg{(j)) in the wire. In Ref. [2l|, a good 
fit for Qe-ph is obtained in the case (/) = by using 
T* — 3.7ETii/kB- The elevated temperature in the wire 
can then be determined by solving 

T!=TU + e^'/^^§^ (6) 

numerically. 

The dominant contribution to the detector noise comes 
from the thermal fluctuations between the electron and 
phonon systems in the N bridge and we refer to this as 
the thermal fluctuation noise (TFN). Due to the super- 
conducting gap, energy exchange between the N bridge 
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and the leads becomes significant only at temperatures 
ksTe ~ A. We first give a simple way to calculate the 
NEP of the bolometer in linear response, i.e., with re- 
spect to small deviations between Te and Tbath- This 
is the usual procedure in determining the noise, valid 
when the input power is small and the system is close 
to equilibrium so that Te w Tbath- In equilibrium, the 
fluctuation-dissipation theorem states that NEPtfn (di- 
rectly related to the thermal noise correlator Sq) is given 
by 



NEP- 



TFN 



therm 



NEpNp^e-^'/^^^ 



Here the thermal conductance of the heat link is 



(7) 



G 



therm 



age 



-ph 



and the effect of superconducting correlations is seen as 
reduction of noise from its normal state value NEP^pj^ by 
an exponential factor. We have also performed a numer- 
ical simulation to determine Gthorm exactly and in Fig. [6] 
we compare this to the two possible exponential formulas 
that can a priori be anticipated to determine the suppres- 
sion of thermal conductance. We see that while the value 
T* = 3.7i?Th/^B used to fit Qe-ph in Ref. [2l|is adequate 
for our purposes in the range ii^Th < fcsTe < 20i?Th, 
the fit is not optimal. Therefore we have also tried 
making the well-founded assumption that the suppres- 
sion of Gthcrm is determined by the niinigap, so that 
Gthorm ^ e:icp{-Eg/kBTe) = exp (-S.lE^Th/fcsre) (at 
(j) — 0). In this case, the fit is improved at temperatures 

Depending on the input optical power, the PJS can 
have Te markedly elevated from Tbath and the linear- 
response result given above is no more sufficient. In this 
case the TEN amplitude has to be calculated from the 
full analysis of the out-of-equilibrium response. We have 
performed this calculation for the noise, yielding 



de de'- 



(e-6'rsgn(e-e') 



24C(5)' 

xK{e,e'){n{e-eV{e'){l-f{e)) 
+ (l + n(e-e'))/(e)(l-/(0)}, 



(8) 



where the Riemann zeta function C(5) ~ 1.0369 
and K(e, e') describes the effect of superconductivity, 
see Ref. |2ll Equation ^ is valid in full nonequilibrium 
where the electron and phonon distribution functions 
/(e) and n(e) can be arbitrary, as long as the electron- 
hole balance is satisfied, i.e., /(e) = 1 — /(— e). We will 
study the nonequilibrium corrections to these results in 
the future whereas here we use the hot-electron model 
where the electron and phonon systems can each be de- 
scribed by a thermal distribution with temperatures Te 
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FIG. 6: (Color online) Comparison between numerical result 
for thermal conductance Gthorm of the electron-phonon heat 
link and two exponential fitting curves: one valid to good 
degree of accuracy for Q (green) and one, where the suppres- 
sion of Gthcrm is determined by the energy minigap Eg (red). 
Gthcrm is given in proportion to the thermal conductance Gjv 
in the incoherent case. 



and Tbath- The distribution functions are given by the 
Fermi function /(e) = (cxp(e/fcBTe) + 1)~-^ and the Bose 
fimction n(e) = (exp(e/A:BTbath) — 1)~^ for electrons and 
phonons, respectively. The structure of Eq. ([S|) can be 
understood by looking at the processes visualized on the 
second and third lines: they correspond to absorption 
and emission of a phonon, respectively. For noise, these 
processes, where heat flows either into or out of the N is- 
land, are additive and consequently show up as a simple 
sum. 

We first estimate analytically the high temperature 
behaviour of Sq from Eq. ([5]). This is done with the 
absorber in the incoherent limit, when the supercon- 
ducting corrections to the absorber are neglected and 
K(e,e') = 1. Equation ^ simplifies to 



305) 



ksT^ J dx x^sgn: 
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Tbath 



1 



(9) 



and, at the two opposite limits of low and high electron 
temperature, this yields: 
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e J 



Te w Tv 



bath 



Sq = ^Mio^nkBTl 



9.8^nkBT' 
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Te 2> Tbath) 
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where the first Hne is in accordance with the linear- 
response result Eq. ([7]). Moreover, both of these limiting 
values are in accordance with the results of Ref. [13 for 
TFN in a nonequilibrium setting with normal conductor. 
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FIG. 7: (Color online) Nonequilibrium noise of the proxim- 
ity Josephson sensor under heating (fcsTbath = 2i?Th), with 
(coherent - red) and without (incoherent - blue) the correc- 
tions from the proximity effect. We compare the exact re- 
sults to the analytical approximations of Eqs. (|10p (incoher- 
ent) and Eq. (O (coherent) valid at the limits Te = Tbath and 
Te ^ Tbath (dashed lines). In the inset, we show how the 
noise of the coherent system can be approximated by the ex- 
ponential law of Eq. ([7|) (dashed line) for a large temperature 
range. 

We also solve Eq. (|8|) numerically to analyze the de- 
tector performance in the coherent limit, where the su- 
perconducting corrections are taken into account. The 
nonlinear noise described by Eq. (|8]) is shown in Fig. [7] 
for our proximity SNS detector compared to the result 
in the absence of the superconducting proximity effect. 
Note that, due to the strong temperature dependence 
of S'q, the limit of high temperatures in Eq. (llObp is 
achieved easily. Already at = 2rbath, the difference 
between the estimate of Eq. ()10b| and the exact value 
from Eq. © is less than 2 percent. This means that 
we can safely use Eq. (jlObp as the upper limit for noise 
whenever the absorber is heated. In addition, as expected 
from Eq. (O, we observe reduction in the noise due 

to superconducting correlations. Nevertheless, in the op- 
erating range of the detector ksTe > -Bxh this reduction 
is relatively small, namely a little more than a factor of 
2 for our given values ksTe = fcsTbath = 2iJTh as seen 
in the inset of Fig. [71 From the inset, we also see that 
the numerically calculated noise follows the exponential 
law of Eq. ([7]) with satisfactory accuracy. This, in turn, 
justifies the use of approximate expression introduced in 
Eq. ([T]) in our analysis. Below, we use a combination 



of Eqs. ([7]) and ([8]) to give analytical estimates of the 
detector performance. 




r(K) 

FIG. 8: (Color online) Base NEP (with negligible power 
load) for proximity and normal conductor as function of 
T = Te — Tbath- The proximity effect decreases NEP by 
the factor ve^^^^Th/s^sT ~ 4 in our example at T = 100 mK. 

Based on the limits of Eqs. (|10p . we may also esti- 
mate the behavior of the detector at linear response, 
i.e., for small input power, and the corresponding NEP 
value for different Tbath is shown in Fig. HI This is the 
base NEP. As seen from Eqs. ([6]) and the heating 
of the absorber due to the incoming optical power starts 
to make a comparable contribution to the noise when 
Pt - ErJT^ji^e^*/^"-"- w 10-18 W... 10-12 W, in the 
temperature range 0.1 K ... 1 K. Beyond this level of 
incident power, the noise becomes independent of Tbath 
and 

NEP ^ yiofc^ \7 ^ « 2 X 10-1° X P^/^ (11) 

in dimensions of W/VHz. For an input power 5 pW, typ- 
ical to balloon and ground-based astronomical telescopes, 
this translates to NEP ~ 3 x 10-^'^ W/VBzM We stress 
that for this type hot-electron bolometer, the noise level 
is almost independent of any other parameters than the 
input power. A factor of 1000 change in the volume of 
the detector changes the noise only by a factor of 2. Only 
when the input power is below ~ TtflT^^^^, which corre- 
sponds to the intrinsic thermal fluctuations of the island, 
can the overall noise be reduced by decreasing the size 
of the absorber. When Pt is larger so that dominant 
contribution to noise results from the radiation-induced 
heating, the noise can actually be reduced, albeit only 
slightly, by increasing the size of the island. 

From the above considerations, we determine the dy- 
namic range of the detector, i.e., the range of input pow- 
ers that can be detected with the sensor. The lower limit 
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is given by NEP and the upper limit, proportional to 
NEP, results from the large increase in temperature Tg 
when the sensor is subject to high level of incident radi- 
ation. The results of increasing Tg are two-fold: firstly, 
the critical current is decreased below the level that can 
be detected by the readout circuit and secondly, the re- 
sponse dIc/dTe becomes too small. The change in tem- 
perature is estimated by taking incident power equal to 
the NEP given in Eq. ^ multiplied by a factor a, and 
using the result of Eq. © at large Tg so that 



and /(e) = [1 -f exp (e/fc^Te)] . Without proximity 
effect, the density of states is independent of energy, 
N{t) = Nf, and'cf = 2t:'^ Np^k^Te/i. The proxim- 
ity effect modifies the density of states by opening the 
minigap Eg ~ iJxh ^ Ai^^ As a result, the heat capaci- 
tance gets modified, but only sligthly, as shown in Fig. |9l 
The growth in Ce is exponential when ksTf, < i?Thj af- 
ter which S saturates and Ce becomes linear, close to the 
normal state value. 



vNEP(Te = Tbath)\/57 
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bath 
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where Sf is the frequency band used in the measure- 
ments. Again, using the values given in Sec. IIIBI and 
assuming Sf = 50 Hz, Tb = 1 K, a = 10^, we obtain 
Te ~ 2.5Tbath, which can realistically be measured. If 
we reduce temperature to Tbath = 0.1 K, w 6.3Tbath 
and for this, one needs to operate close to Tbath = T-Th so 
that Te does not become too large. As the dynamic range 
is considered, it remains important to decide whether we 
wish to obtain linear output from the detector, so that 
only a restricted temperature range is available due to 
the exponential depence of Ic at high temperatures (as 
in Fig. The dynamic range is therefore dependent 
on how the detector is optimized but typical values of 
a = 10^ are still well within reach. 



Calorimeter 




FIG. 9: (Color online) Comparison between the numerical re- 
sult for the electronic heat capacity Ce of the coherent (prox- 
imity) and the incoherent case. 



In calorimetric mode, the rate of incident photons is 
much smaller than the rate of recovery for the detector. 
Consequently, photons are absorbed one at a time and 
after each absorption event, Te rises rapidly within the 
diffusion time td = D/P, and then decays back to the 
equilibrium value Tbath within time r, determined by the 
relaxation processes. We assume that the rise time of the 
pulse is very short, tjj <^ r, compared to the thermal 
time constant of the detector, which reads 



Ce 



Gth crm 

The heat capacity of the electron system is given by 

dS{Te) 



(12) 



CeiTe)=Te 



dTe 



(13) 



where the entrop; 



,24.25 



S ^-2kBfl J deN{e) 



(14) 



X {/(e)ln [/(6)] + [1-/(6)] In [1- /(e)]}.. 



Now that we have solved for the heat capacity C and 
thermal conductance Gthcrm numerically, the time con- 
stant T may be obtained from them using Eq. (jl2p . Since 
the heat capacity as seen in Fig. [9] is approximately lin- 
ear, we may conclude that the temperature dependence 



of T follows the law t 



^*/T/y3^ The exact resuhs 
from numerical simulations are shown in Fig. [TU] and from 
there, we see that rc_ph ^ 10^"^ . . . 10~^ s in the range 

Te = 0.1...1.0K. 



The calorimetric mode is characterized by the resolving 
power 



hv 

Je 



2V21n 2NEPtfnv^ 



The estimate for the background noise, Ai?, should be 
based on the idle state of the detector, i.e., on the state 
where Tg = Tbath and no signal is present. Here we rectify 
the erroneous statement made in Ref. 3, where the noise 
is overestimated and the resolving power correspondingly 
underestimated. We may use the approximative formulas 
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FIG. 10: (Color online) Time constant of the proximity 
Josephson sensor as a function of electron temperature on 
a log-log scale. The coherent corrections from the proximity 
effect (red) are compared to the incoherent result t ~ T^"^ 
(blue). 



given above for NEPtfn and r to obtain 

hv hv 

SE ~ 4v/21n2fci3T4ti^GthcrmV^ 

_ hv j 3 

3/2 

The scaling with T^g^^{^ is approximatively valid for the 
temperature scales of our interest, where Ce is nearly 
linear. Energy resolution of unity is obtained when 

= 8n^ln2Npnk%TU3h^ « 16 THz x f^i^^, (16) 

and Tbath is the bath temperature in units of Kelvin. This 
means that at 100 mK, it should be possible to realize a 
single-photon detector for frequencies above 500 GHz. At 
50 mK, this limit is 180 GHz. In contrast to the bolome- 
ter case, there is some freedom to alter the numerical 
prefactor in Eq. (jl6l) by choosing different material pa- 
rameters from those given in the beginning of Sec. IIIII 
Since the proximity effect modifications in r and NEP 
are in opposite directions, the resolving power is almost 
the same as that given in Eq. (IT5t . which neglects the 
proximity corrections. This is shown in Fig. 1111 



FIG. 11: (Color online) Detector resolving power for some 
operating temperatures. Note that the curves for Tbath = 
1.0 K are nearly on top of each other. 



C. Optimization 

The scheme described above is largely adjustable and 
such adjustment is required whenever the objective of the 
detector operation is changed, for example in the case of 
different power loads Pt or operating temperatures Tbath- 
While the exact optimal parameters depend largely on 
the purpose and the situation where the detector is used, 
we now provide some of the most interesting examples 
of the effects of altering the exact detector configuration 
and materials. 



1. Length 

The most notable material-independent characteristic 
of the detector is the absorber length This affects the 
Thouless energy E'xh ^ and, consequently, also the 
relation between Ic and Te, as shown in Fig. [5] and the 
corresponding Eq. ([T]). Since the detector is based on 
the strong observable response in Ic, which again is due 
to changes in Te, altering Exh markedly shifts the elec- 
tron temperature regime where the detector operates. In 
Sec. nil Al it is shown that in bolometric operation it is 
either the intrinsic temperature fluctuations or the inci- 
dent power itself that governs the noise level. Therefore 
it becomes crucial to be able to find the most appropriate 
operating temperature regime. 

The operating regime is determined for two parameters 
Te and Tbath so that T^in < ?bath < ^e < Tmax for the 
lower and the upper limits of the regime, Tmin and Tmax- 
In Sec. Ill 21 we determine the lower limit of operation 
to be fcsTmin = l-Oi?Th sincc below this temperature, Ic 
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saturates. The upper limit is set by our approximation 
of the minimum measurable critical current 3 nA. For 
our base value of absorber length, I = 1 fim, this cor- 
responds to T,„ax = 20i?Th — 400 mK. The change in 
length has two effects. When i?Th becomes smaller, the 
curve in Fig. [2] shifts left and the response diminishes, 
and vice versa. On one hand, long absorbers thus imply 
reduced supercurrent but on the other hand, this also 
means that Erh is reduced and the detector can be oper- 
ated at lower temperatures due to increased response in 
this regime. Conversely, if the absorber is short, the elec- 
tron operating regime of the detector is shifted to higher 
temperatures and the critical current is increased for an 
easier readout. 

When / is increased by a factor of 2, fcsTmin = 
^'Th — >■ jEtyi ~ 5 mK. Due to the decreasing super- 
current, we now have Tmax ~ 55 niK. These tempera- 
tures are already too small for present-day applicability. 
On the other hand, when I is decreased by a factor of 2, 
ksTmin = Eti, ^ 4Ety. w 76 mK and T^ax ~ 2.4 K. 
The drawback of this large temperature range is that 
we can expect the outdiffusion of energy to increase for 
the short absorber. This would change the fraction of 
the radiation power r{i') which yields a thermal signal 
as discussed in Sec. Ill Al We assess how changing of the 
absorber length affects r{iy) by noting that the only term 
in Eq. ^ which scales with I is the diffusion part pro- 
portional to ^ Changing / therefore only changes 
the relative strength of diffusion compared to the irradi- 
ation and inelastic relaxation. The contribution to r{uj) 
due to diffusion comes from energies close to the gap A, 
as seen in Fig. |31 However, in our parameter regime the 
effect of diffusion in Eq. ^ is so weak compared to the 
inelastic energy re-distribution I^-c and /c-ph that a fac- 
tor of 2 change in length of the absorber results only in 
a change of the order of 10-20 percent in r at low fre- 
quencies hiy < 3A. At large frequencies, the energy re- 
distribution through inelastic collisions is the dominant 
mechanism of power loss. We therefore conclude that the 
amount of radiation power that is coupled to the sensor 
is relatively independent of small changes in / in our pa- 
rameter regime. 



2. Graphene 

In Sec. nil A[ we have also determined that reducing 
the volume of the absorbing element, the intrinsic ther- 
mal noise can be decreased (see Eq. ([7])). Consequently, 
graphene, which is composed of a single layer of graphite, 
can be seen as an attractive alternative for the absorber 
material when extremely weak radiation signals such as 
single photons are detected and Popt is small. It has al- 
ready been demonstrated that graphene, placed between 
superconductors, can support large supercurrenta^ii^ 
For graphene, the heat flow between the electrons and 
phonons does not follow Eq. ^ anymore since this for- 
mula is valid only for 3D conductors. In Ref. [2^, the 



electron-phonon cooling power in graphene has been de- 
termined at temperatures with a minimum of = 20 K. 
While this is already far above our operating regime, we 
may calculate Qo-ph for our metallic wires at such tem- 
peratures to get some idea where graphene stands when 
these two different realizations are compared. 

According to Ref. [2^ a graphene flake of width w = 
9 /im and length I = 350 nm with gate voltage such 
that carrier density becomes n ~ 10^^ . . . 10"'^'^ 1/cm^ has 
prominent characters such as diffusion coefficient com- 
parable to the values of our example above. More im- 
portantly, for such gating it has normal-state resistance 
Rn ^ 50 O and could thus be used potentially as an ab- 
sorber with good impedance matching. Since the length 
of the graphene absorber is reduced to one third of the 
value used for our metallic absorber, i?Th ^ is cor- 
respondingly increased by a factor of 9 and the operat- 
ing regime in temperature shifts from 20 mK. . . 400 mK 
to 200 mK. . . 3.5 K. We use these specifications for the 
graphene detector. Then, from Ref. [2^ the electron- 
phonon cooling power becomes 300 nW. . . 3 fiW [n ~ 
10"^^ . . . 10"^^ 1/cm^). This should be compared with the 
value we get for the electron-phonon cooling power when 
we use Eq. (|4]) for our metallic conductor at Te = 20 K 
and Tbath = 0. By extrapolating outside the operating 
regime of the metallic detector this way, Qo-ph ~ 1.6 /zW. 
This is in the middle of the values we have for graphene, 
and we conclude that graphene can be used as a tun- 
able detecting element with gate voltage determining the 
strength of relaxation. It is difficult to estimate what 
happens when the temperature is around = 1 K since 
direct extrapolation from the results of Ref. |2^ would 
require knowing the exact temperature dependence of 
Qc-ph for graphene. However, we expect that by adjust- 
ing the carrier density in graphene by the gate voltage, 
the thermal conductance of the electron-phonon link can 
be made small in comparison to the case of the metallic 
wire, thus producing favourable NEP figures. 



IV. CONCLUSION 

In this paper, we present in detail the operating con- 
cept of the proximity Josephson sensor introduced in 
Ref. and study the physical processes involved in the 
detector operation. First, we devise a model to explain 
how the incident radiation couples to the sensor, how the 
radiation energy is distributed inside the absorber after 
the initial excitation and how the energy finally escapes 
the system. We adopt the hot-electron model where the 
absorption of radiation results in increase of the electron 
temperature. Then we show how the superconducting 
proximity effect in the absorbing element of the detector 
alters the heat capacity and the heat conductance of the 
detector and how this affects the detector performance 
in terms of the nonequilibrium noise. We observe that 
in the operating regime of the detector, where temper- 
ature is larger than the Thouless energy, corrections to 



12 



the expected magnitude of noise due to the proximity 
effect are small. Nevertheless, they reduce the intrinsic 
thermal noise of the detector, which sets the limit for 
detector performance. 

Since the detector is largely adjustable depending on 
the operating temperature and power load, we refrain 
from giving any definite limits of performance. How- 
ever, we show the physical limit of minimum noise for 
a given set of material parameters in terms of the base 
NEP in Fig. [51 This base NEP sets the ultimate limit 
for single-photon detection in the calorimetric mode and 
in the bolometric mode it sets the limit for small input 
powers. For large input powers, this NEP is to be com- 
pared with the NEP arising from the heating, given in 
Eq. ([Tl]) . Notably, the latter is only weakly dependent 
of the volume or the strength of electron-phonon cou- 
pling in the absorber. We observe that at T ~ 100 niK, 
base NEP of the order of 10^^° W/\/Hz can be expected. 
Under power loads of several picowatts, corresponding 
to a temperature raise to a few Kelvin, NEP becomes 
~ 10"-'^^ W/-\/Hz. For power loads this large it is de- 
sirable to fabricate devices with larger volumes, leading 
to smaller temperature changes and, consequently, lower 
NEP. It follows that when designing a detector for a par- 
ticular purpose using the PJS concept, we need to bal- 
ance between the intrinsic (base) noise of the absorber at 
Tf, — Tbath and the noise that results from heating the 
absorber above Tbath- 

The dynamic range of the detector can be of the order 
of 10^ but if the detector is operated at a large range of 
input powers, the nonlinearity in response needs to be 
taken into account. On the other hand, the numbers for 
detector resolving power suggest that our scheme could 



also be used effectively in detection of extremely weak ra- 
diation signals. This includes single-photon detection for 
frequencies above 180 GHz and even detailed radiation 
spectroscopy at the terahertz regime. Optimizing the 
detector for a specific task requires choosing appropriate 
materials, and for this we have studied the behavior of 
graphene as the absorber. Notably, both the graphene 
and metallic wire concepts share the advantage that, by 
controlling the length-to-width ratio of the absorber, a 
good impedance matching to the antenna circuit is pos- 
sible. However, we encounter additional mechanisms of 
power loss in optical coupling of the incident power to the 
absorber, namely the outdiffusion of quasiparticles above 
the superconducting energy gap and nonthermal redis- 
tribution of radiation-excited quasiparticles. In a typical 
range of materials parameters, the latter is so dominant 
that the length of the absorber may be varied without 
experiencing a drastic increase in outdiffusion allowing 
the detector to be optimized for varying operating re- 
quirements. This is especially important if the detector 
is to be operated at temperatures over 1 Kelvin. In that 
case, it needs to be investigated whether sub-micron or 
graphene absorbers may be used. 
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